#!/usr/bin/python
# This script will parse and print Fasta files >= a given size

from Bio import SeqIO
from Bio.SeqRecord import SeqRecord
import sys

handle = open(sys.argv[1], "rU")

i = str(sys.argv[2])

records = list(SeqIO.parse(handle, "fasta"))

for rec in records:
    iden = rec.id
    seq = rec.seq
    if iden == i:
        print ">" + iden
        print seq

handle.close()
